Renormalization and cyclotron resonance 
in bilayer graphene with weak electron-hole asymmetry 
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Cyclotron resonance in bilayer graphene is studied with weak electron-hole asymmetry, suggested 
by experiment, taken into account and with the focus on many-body corrections that evade Kohn's 
theorem. It is shown by direct calculation that the theory remains renormalizable to 0{e'^) in the 
presence of electron-hole asymmetry parameters, and a general program to carry out renormalization 
for graphene under a magnetic field is presented. Inclusion of electron-hole asymmetry in part 
improves the theoretical fit to the existing data and the data appear to indicate the running of the 
renormalized velocity factor with the magnetic field, which is a key consequence of renormalization. 

PACS numbers: 73.22.Pr,73.43.Lp,76.40.+b 



I. INTRODUCTION 

Graphenei^— supports charge carriers that behave as 
Dirac fermions, which, in a magnetic field, lead to a 
particle-hole symmetric and unequally-spaced pattern of 
Landau levels. Accordingly, graphene gives rise to a 
variety of cyclotron resonance, both intraband and in- 
terband, with resonance energies varying from one res- 
onance to another. This is in sharp contrast to stan- 
dard quantum Hall systems with a parabolic energy 
dispersion, where cyclotron resonance takes place be- 
tween adjacent Landau levels, hence at a single frequency 
(jjc = eB/m*, which, according to Kohn's theorem^- is 
unaffected by electron-electron interactions. The non- 
parabolic spectral in graphene offer the challenge of de- 
tecting many-body corrections to cyclotron resonance. 

Theoretical studies^"— over the past few years have 
revealed some notable features of quantum corrections 
to cyclotron resonance in graphene and bilayer^i^ill 
graphene. The genuine many-body corrections arise 
from vacuum polarization, specific to graphene, which di- 
verges logarithmically at short wavelengths. This means 
that one has to carry out renormalization properly, as 
in quantum electrodynamics, to extract observable re- 
sults. In particular, for bilayer graphene it turns out 
that both the leading intralayer and interlayer coupling 
strengths undergo renormalization and that their renor- 
malized strengths run with the magnetic field. Bilayer 
graphene is marked with the unique property that its 
band gap is externally controllable. ^^"^^ 

Experiment has so far verified, via infrared spec- 
troscopy, some basic features of cyclotron resonance in 
monolayer—"— and bilayeri22,i2l graphene. The data for 
the monolayer show a good symmetry between the elec- 
tron and hole bands but generally show no clear sign of 
the many-body effect, except for a datumf22, Indeed, a 
comparison between some leading intraband and inter- 
band cyclotron resonances revealed a small deviation in 
excitation energy, consistent with the presence of many- 
body corrections roughly in magnitude and sign. 

The situation is quite different for bilayer graphene. 



for which only a limited number of data are available so 
far. The dataS^ on intraband resonances show a weak 
electron-hole asymmetry, and generally defy a good fit 
by theory. Actually one has to employ different values 
of the velocity factor v to fit the electron data and hole 
data separately. 

Earlier Raman^^ spectroscopy and subsequent in- 
frared^^"— spectroscopy of bilayer graphene under zero 
magnetic field also revealed a significant asymmetry be- 
tween the conduction and valence bands, mainly due to 
subleading intra- and inter-layer couplings A and 74. 

The purpose of this paper is to reexamine cyclotron 
resonance in bilayer graphene, with possible electron-hole 
and valley asymmetries taken into account. It is not clear 
a priori whether the renormalizability of the low-energy 
effective theory is maintained in the presence of electron- 
hole asymmetry, since the asymmetry parameters (espe- 
cially, 74) critically modify the ultraviolet structure of the 
theory. We show that the theory indeed remains renor- 
malizable to O(e^) (at least), and that the renormaliza- 
tion counterterms depend on 74 in a nontrivial way. We 
present a general algorithm to carry out renormalization 
for graphene under a magnetic field, executable even nu- 
merically. Inclusion of electron-hole asymmetry parame- 
ters partially improves the theoretical fit to the existing 
data, and the fit in turn suggests some nontrivial modifi- 
cation of the spectra of the zero-mode and pseudo-zero- 
mode Landau levels specific to bilayer graphene. 

In Sec. II we briefly review the effective theory of bi- 
layer graphene and examine the effect of electron-hole 
and valley asymmetries. In Sec. Ill we study the Coulom- 
bic many-body corrections to cyclotron resonance, with 
a focus on renormalization and its consequences. Sec- 
tion IV is devoted to a summary and discussion. 



II. BILAYER GRAPHENE 

Bilayer graphene consists of two coupled honeycomb 
lattices of carbon atoms, arranged in Bernal A'B stack- 
ing. The electrons in it are described by four-component 
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spinor fields on the four inequivalent sites {A, B) and 
{A',B') in the bottom and top layers, and their low- 
energy features are governed by the two inequivalent 
Fermi points K and K' in the Brillouin zone. The in- 
tralayer coupling 70 = jab ^ 3 eV is related to the 
Fermi velocity v = (\/3/2) aL7o/?i ~ 10^ m/s (with al = 
0.246nm) in monolayer graphene. The interlayer cou- 
plings2^i^ 71 = jA'B ~ 0.4 eV and 73 = jab' 0.1 eV 
are one-order of magnitude weaker than 79. Actually, 
interlayer hopping via the {A',B) dimer bonds modifies 
the intralayer linear spectra to yield quasi-parabolic spec- 
tra^** in the low-energy branches |e| < 71. 

The bilayer Hamiltonian with the leading intra- and 
inter-layer couplings w cx 70 and 71 lead to electron-hole 
symmetric spectra. Infrared spectroscopy^ of bilayer 
graphene, however, has detected some weak asymmetry 
between the electron and hole bands, such as (i) the en- 
ergy difference A « 18 meV between the A and B sub- 
lattices within the same layer and (ii) the next-nearest- 
neighbor interlayer coupling 74 = jaa' — Ibb' ~ 0.04 70. 

The effective Hamiltonian with such intra- and inter- 
layer couplings is written as^^'^^ 



hi 



dh 



V p 



Vip vp 



t A- 

\ vp V4P^ 71 



up' 

V4P 

71 



A 



i„ 



(2.1) 



with p — px + ipy, — Px — iPy Here ^I' = 
(^A, V'-B' J V'A' 7 "05)' stands for the electron field at the 
K valley, with A and B referring to the associated sub- 
lattices; u stands for the interlayer bias, which opens a 
tunable gapiS, between the K and K' valleys. We ig- 
nore the effect of trigonal warping cx 73 which, in a 
strong magnetic field, causes only a negligibly small level 
shift. ■^^ We also ignore weak Zeeman coupling and, for 
conciseness, suppress the electron spin. Our definition of 
V4 = —(74/70) V differs in sign and by factor v from the 
one (^4 — >■ 74/70) in the literature^^""— i^^; this choice is 
made simply for notational convenience. 

The Hamiltonian H- at another {K') valley is given by 
'H+ with {v, V4, u) — ?> (—V, — W4, —u), and acts on a spinor 
of the form ^ = {^pB' ,'4'A,tpB,tpA'Y ■ Note that 7^+ is 
unitarily equivalent to "H- with the sign of u reversed. 



(2.2) 



with U = diag(l, 1,-1,-1). This implies that the elec- 
tronic spectrum at the K' valley is obtained from the 
spectrum at the K valley by reversing the sign of u; in 
particular, the spectra at the two valleys are the same 
for M = 0. Nonzero interlayer voltage u ^ thus acts as 
a valley-symmetry breaking. 

We adopt the set of experimental values-- 

u sa 1.1 X 10^ m/s, 71 w 404 meV, 
V4/V = -74/70^-0.04, A«18meV, (2.3) 



in what follows. Full account is also taken of the effect 
of interlayer bias u. For notational simplicity, however, 
we often present analytical expressions only for u — 0. 

The Hamiltonian H^' gives rise to four bands with elec- 
tron and hole spectra, which, for u — 0, read 



e4(p) = 

e3(p) = 

e2(p) = 

ei(p) = 



«|p2 + (7+/2)2+7+/2. 



p2 + (7_/2)2-7_/2, 

V^'iP' + (7+/2)^+7+/2, 

p2 + (7_/2)2-7_/2, (2.4) 



where v± = v ± V4 and j± = 71 ± A. Note that V4 and 
A effectively modify v and 71 , respectively, in a manner 
different for electrons and holes; the spectra are electron- 
hole asymmetric unless D4 = A = 0. These band spectra 
acquire nonzero valley gaps for u 7^ 0. 

Let us place bilayer graphene in a strong uniform mag- 
netic field Bz = —B < normal to the sample plane; 
we set, in H±, p ^ U = p + eA and p'' — >■ H''' with 
A = Ax+iAy — By, and denote the the magnetic length 
as ^ = 1/^/eB. It is easily seen that the eigenmodes of 
'H+ have the structure 



n-l)bif\\n-l)bit^ 



(2.5) 



with n = 0, 1, 2, . . ., where only the orbital eigenmodes 
are shown using the standard harmonic-oscillator basis 
{\n)} (with the understanding that |n) = for 7i < 0). 

The coefficients b„ = {b'i^\ b'k\ b'i\ bl^^ f for n = 2, 3, . . . 
are given by the eigenvectors of the reduced Hamiltonian 



-In' 



7' 

d+\u' j 



(2.6) 



where 



\/2v/£ « 36.3 X t>[10V/s] ^^[T] meV, (2.7) 



with V measured in units of lO^m/s and B in tcsla, is the 
characteristic cyclotron energy for monolayer graphene; 
r = V4/V = -74/70 (~ -0.04), 7' = 7i/wc, d = A/wc 
and u' = u/ujc- 

The energy eigenvalues e„ of "Hrod are determined from 
the secular equation, which, for m = 0, reads 

n{n - 1)(1 - r^ f - ~ d' + c„ (1+ r')] e'' + {e'^f 



-c„[2r7'-d(l+r2)]e'-2fi {(') 



/n3 



0, 



(2., 



with c„ = 2?i — 1 and e' = tnjujc- We first consider the 
M = case. Let us denote the four solutions of the secu- 
lar equation as el„ < e_„ < < e„ < e+ for each integer 
n > 2, so that the index ±n reflects the sign of the en- 
ergy eigenvalues. For n = i/^' has an obvious zero 
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eigenvalue eo = 0, with the eigenvector bo = (1,0,0,0)*. 
For n = 1 the secular equation (|2.8p is reduced to a 
cubic equation in e', excluding e' = 0, and leads to 
three solutions, which, for the present choice (|2.3p of pa- 
rameters, are given, e.g., by e' — (—3.96,0.029,4.29) at 
_B = 10 T; we thus denote the corresponding eigenval- 
ues as (el]^, ei, £]*"). This eigenvalue ei changes sign if 
one sets d — —d and r — > — r simultaneously. Thus the 
assignment of e±i in general depends on the choice of 
asymmetry parameters (v4, A) and also on the interlayer 
bias u. Actually, for zero bias m = 0, ei deviates from 
zero as V4 and A develop. In this sense, the n = 1 Lan- 
dau level is a pseudo-zero-mode level while the n = 
level is a genuine zero-mode level. 

As u is turned on, these n = (0,±1) levels go up or 
down oppositely at the two valleys; e.g., cqIk = —£o\k' = 
u/2. Their spectra vary linearly with u while other lev- 
els (|n| > 2) get shifted only slightly. Interestingly, for 
< |u| ^ Wc, ei ^ eo^. at one valley while e\n\=i at 
another valley crosses^E eg^ from below with increasing 
magnetic field B. A nonzero u thus critically spoils the 
valley symmetry of the n = (0, ±1) sector. 

The spectra and el„ with n > 1 form the 
high-energy branches of the electron and hole Lan- 
dau levels, respectively; |e^| > 71. Let us com- 
bine £±n into the low-energy branch of Landau levels 
{cn} = {••■! £-3: £-2, eo, ei, 62, ■ • ■}, and denote the three 
branches (e~ 



en ^ cjcVn{l',d,r,u'), ^ oJcVnil' id, r,u'). (2.9) 

These {rjn.'nn) ~ ^' are uniquely determined from 7^ red 
of Eq. (|2.6p or from Eq. (12. 8p as functions of (7', d, r,u'). 
One can thereby construct the associated eigenvectors, 
which, e.g., for u = and \n\ > 1, read b„ = 



{b. 



(1) ,,(2) ,(3) ,(4).t _ ^(1) n(2) 



with 



«(3) 
h'n 

«(4) 
h'n 

«(2) 
h'n 



\n\ - Tin (?7„ -d) + (|n| - 1) r'^ 



n\ Gn 

7'?7„ + (2|n|-l)r]/(v^G„), 
/H^(/3i'^+r/3i'^ (2.10) 



where G„ = 7'+r {7?„-d+(|n|-l) (l-r^)/77„} and 6^^^ = 
1/^ 



=2,3,4 



(On 2 



These expressions are equally 



valid for eigenvectors belonging to e„ , with rjn — ^ r]i^ . 

The Landau-level spectrum e„ — uJcVniy ,d,r,u') de- 
pends on the magnetic field i? in a nontrivial manner 
through the dimensionless quantities 7', d and u'. Ac- 
tually, for the choice of 71 = 404 meV, A = 18 meV, 
r = —0.04 and u = 0, the electron and hole spectra differ 
considerably, as shown in Fig. 1. In particular, the zero- 
mode level {n — 0) remains intact while the pseudo-zero- 
mode level {n = 1) gets shifted, e.g., by ~ 8 meV as B is 
increased from to 20 T. The Landau gaps are generally 
larger for electrons than holes; this is readily understood 
from the behavior of e3(p) w (w^/7'_)p^ for p — > 0, 




FIG. 1: The 3 < 2 (or 2 < 3) resonance (a) Landau levels 

as a function of magnetic field B for 71 = 404 meV, A = 18 
meV, V4/V = —0.04 and = 0; e„ for electrons (solid curves) 
and —e-n for holes (dashed curves), (b) Landau gaps for 
electrons (solid curves) and holes (dashed curves), (c) Typical 
channels of cyclotron resonance for u = 0; circularly polarized 
light can distinguish between the two classes of transitions 
indicated by different types of arrows. 
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FIG. 2: Landau gaps in Fig. 1 (b), compared with the ex- 
perimental data in Ref. [S^ on cyclotron resonance energies, 
(a) Electron band, (b) Hole band. 



which implies that the effective mass m* ~ 7_/(2t)^) is 
smaller for electrons. The asymmetry between the elec- 
tron Landau levels L„ and hole levels L_„ becomes more 
prominent for higher levels |n| > 2. 

Cyclotron resonance in bilayer graphene is governed by 
the selection rule^ A|n| = ±1, and there are two classes 
of transitions, (i) i±(n-i) -^-n and (ii) L„ -s— -/j±(n_i) 
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(with n > 2), which are distinguished by the use of 
circularly-polarized light. See Fig. 1 (c). 



III. CYCLOTRON RESONANCE AND 
MANY-BODY CORRECTIONS 



As for experiment, Henriksen et al^ measured, via 
infrared spectroscopy, cyclotron resonance in bilayer 
graphene in magnetic fields up to 18T. They observed in- 
traband resonances, which are identified with the L2 ^ 
Li, L3 ^ L2, L4 L3 and L5 ^ L4 transitions at 
filling factor i'— 4, 8, 12 and 16, respectively, and the 
corresponding hole resonances at 1^= -4, -8, -12 and -16, 
together with a significant asymmetry between the elec- 
tron and hole data. 

Such data from Ref. [l^ are reproduced in Fig. 2. Also 
included are the zeroth-order Landau gaps of Fig. 1 (b), 
which apparently fit the experimental data reasonably 
well, except for the v = A electron data on the L2 ^ Li 
resonance, which deviates considerably. 

It is worth discussing the effect of interlayer bias u here. 
The n — (0, 1) levels are very sensitive to u and may eas- 
ily acquire valley gaps for m ^ while other levels |n| > 2 
are relatively inert as long as |m| uJc- Cyclotron res- 
onances involving the n ~ 1 level, i.e., (1 —2) and 
(2 1) resonances, therefore tend to be affected by u. 
Actually, with u ~ 20 meV one can apparently fit the 
data for those resonances at one valley. The asymme- 
try, however, is reversed at another valley. Nonzero u 

may thus broaden the observed widths of the (1 < 2) 

and (2 -s— 1) resonances but would not account for their 
asymmetry. 

It is clear now that one should treat those resonances 
separately from the rest of the resonances, which are 
barely sensitive to |m| <C lUc- The former and latter are 
also different in their sensitivities to electron-hole asym- 
metry a A and V4. See Fig. 1 (b) again. It shows that 
for It = the 2^1 and 1 — 2 gaps barely differ 
while other Landau gaps exhibit significant asymmetry 
between the electron and hole bands. In this sense. Fig. 2 
shows us that A w 18meV and U4 w —O.OAv, obtained 
from independent experiments, account for the electron- 
hole asymmetry between the = (8, 12, 16) data and the 
1/ = (—8, —12,— 16) data reasonably well. It is a non- 
trivial fact that this single set of parameters can fit the 
electron and hole data simultaneously. 

Poor fitting to the f = 4 data, on the other hand, 
would suggest that the spectrum of the pseudo-zero- 
mode level (|n| = 1) is further modified^ by some other 
sources. Actually it is expected theoretically^"— that 
the n — (0, ±1) sector of bilayer graphene has nontriv- 
ial dynamics due to orbital mixing and supports charac- 
teristic collective excitations, orbital pseudospin waves. 
It would be interesting to study how the electron-hole 
asymmetry affects the detailed structure of this special 
sector. 



In this section we study the many-body corrections to 
cyclotron resonance, with emphasis on how to carry out 
renormalization. The Coulomb interaction is written as 



H 



Coul 



2 ^ 
p 



i^P : P-pPp ■; 



(3.1) 



where pp is the Fourier transform of the electron density 
p = 'qj^^ + ^t^; 27ra/(eb|p|) is the Coulomb po- 

tential with a = e^/(47reo) « 1/137 and the substrate 
dielectric constant Cb; J2p — I '^^P/(27r)^. 

The Landau-level structure is made explicit by pass- 
ing to the \n,yo) basis (with yo = i'^px) via the expan- 
sion (*(x,t),v[/(x,t)) = J2n^yo(^\^^yo) '^niyo:t); remem- 
ber that fields carry (suppressed) spin and valley in- 
dices. The Hamiltonian H^^ is thereby rewritten as 

/oo 
dyoY,^iiyo,t)enMyo,t), (3-2) 
n— — oo 

and the charge density p-p{t) = / d^xe^^'^ p as^ 

00 00 

Ekn \ ^ kn inkn 

p-p = 2^ 9p Rp , 

k,n— — 00 k,n— — 00 

^p" = 7p / dyo^Pi{yo,t)e'P-'iJn{yo,t), (3.3) 



where 7p = e ^ p /4- j. — (ii'^d/dyo, yo) stands for the 
center coordinate with uncertainty [rx,ry] = il?. The 



^36 



charge operators R'^ obey two Woo algebrasiiii associ- 
ated with intralevel center-motion and interlevel mixing 
of electrons. 

The coefficient matrix q^;" is constructed from the 
^p 

knowledge of the eigenvectors b„ , 



kn 
Jp 



where 



Jp 



kn 




(3.4) 



(3.5) 



for fc > n > 0, and f^'' = (/:^p)^; p = Px+ipy Ex- 
pression (13. 4p is valid for \n\ = 0, 1 as well, with the 



understanding that /p" = for A; < or n < 0. 

For zero bias u = 0, are the same at the two val- 

leys, i.e., 5^"!^' = 5p"k with {b[l\b':^\b[? ,b^^^)\K' = 

{hn\hn\—hn\—b'"n^)\K- This follows from the unitary 
equivalence p.2p of the Hamiltonians T-L±\±u and the in- 
variance of the charge density pp under U there. 

The Coulombic correction to cyclotron resonance in 
graphene to 0{a/t) was calculated earlier^^ using the 
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single-mode approximation.^^ Here we consider cyclotron 
resonance (at integer filling v) from the filled ath Landau 
level (La) to the empty 6th level (Lf,) at zero momen- 
tum transfer k = 0, where no mixing takes place in spin 
and valley. The cyclotron-resonance energy for a general 
Lb -i— La transition with the Landau levels filled up to 
n = Hf (a < Uf < b) is written as^^ 



= £6 - 



(3.6) 



with the correction 



an 1 2 
pl 



\9p\') 



(3.7) 

diagonal in spin and valley. As shown by Fig. 1 (c), 
Tif — —4,— 3,— 2,1,2,3,... correspond to the filling factor 
v = -12, -8, -4, 4, 8, 12, respectively. 

One can now substitute Eq. p.4[) into this formula 
and calculate the Coulombic corrections with the ef- 
fect of electron-hole asymmetry taken into account. 
There is, however, one technical problem to solve. The 

En<„f (Iffp"!' - Iff-pl") in <E3) refers to quan- 
tum fluctuations of the filled states and actually diverges 
logarithmically with the number Nj^ — >■ cxd of filled Lan- 
dau levels in the valence band (or the Dirac sea). 

One has to handle such ultraviolet (UV) divergences 
by renormalization of the basic parameters (w, W4, 71, A) 
and, if necessary, the tunable parameter u. In the 
electron-hole symmetric case, v and 71 turn out to be 



rcnormalized in the same wayj*^ i.e. 



Zyv"""" and 



71 = 75;™ ^[^^ Zy = to 0(a/i). Actually, it is not 
clear a priori if the theory remains renormalizable in the 
presence of asymmetry parameters U4 and A. If inclu- 
sion of V4 and A were to yield a new type of divergence 
unremovable by rescaling of the existing parameters, the 
theory would lose renormalizability (or one would have 
to introduce a new parameter to remove the divergence 
and, if necessary, repeat this process). We prove by di- 
rect calculations below that the theory is renormalizable 
to 0{a/£). 

The key to this problem of renormalization is to note 
that the magnetic field supplies only a long-wavelength 
cutoff through the magnetic length i — 1/^/eB, leaving 
the UV structure of the theory intact. One can therefore 
first look into the theory in free space {B = 0) and deter- 
mine the UV structure of the Coulomb exchange correc- 
tions. Such corrections are written as fk *>S'(p -I- k), 
a convolution of the photon propagator — 27ra/(ef,|k|) 

FT 

and the instantaneous electron propagator (^'5'' = ' 
iS{p). Their UV structure is thus read from the asymp- 
totic behavior of S{p). 

The resulting divergences are then absorbed into the 
counterterms Sv, Sji, Sv^, SA and Su, generated by 
rescaling 



V 

71 



I, V 



Sv, 



(3.8) 



where "ren" refers to renormalized parameters. See Ap- 
pendix A for such an analysis of divergences. Here we 
quote only the result: (i) The counterterm for velocity 
factor V turns out to be the same as in the electron-hole 
symmetric case (and in the case of monolayer graphene^l 
as well). 



Sv = {Zy - l)v"' 



ia/Seb) logA^ 



(3.9) 



where A stands for the momentum cutoff which is related 
to the Dirac-sea cutoff so that^^ A^ « 2Nl/P- (ii) Re- 
markably, U4 and u remain finite, 



(5w4 = (5u = 0, 



(3.10) 



and require no renormalization, W4 = v''^'^ and u — u^™. 
(iii) The dimensional parameters 71 and A are mixed 
under renormalization, 



= (7r-^'' 

SA = 2(A'™- 



A' 



^rcn ^rcn 



)/i[rn 6v/v' 



(3.11) 



where h[r] = 1/(1 — r^) and 



"/u'-™. Note that 



the counterterms are highly nonlinear in r''™ oc W4. 

One can now pass to the B case with these coun- 
terterms. Let us denote by "H™" the Hamiltonian H+ 
[of Eq. (|2.ip ] in magnetic field B with {v, V4, 71, A, u) re- 
placed by (1;''°", W4™, 7^'^", A''™, u''''"), and write its spec- 

rcn /'^/rcn ^rcn 



trum as e„ 



) with u}l° 



^v'°'^ jl, etc., in obvious notation; see Eq. p.9p . Sup- 
pose now that we start with Hip" and calculate Coulom- 
bic corrections to 0{a/ tj. The divergences we encounter 
are removed by the counterterms formally written as 
(^ct^™! where the differential operator 



5v- 



d 



d 



97- 



5A- 



d 



dA' 



(3.12) 



acts on 'HiJ?". For the related reduced Hamiltonian 'HJ. 
defined as in Eq. (j2.6p , the counterterm is also written as 



en 
d' 



(5ct'HJ.°3- for example, one is to subtract divergences 
from the 0{a/t) correction to the spectrum e™°, the 
required counterterm is obtained from the expectation 
value of Jct^fed' which equals^^ 6ct "> the variation of 
the eigenvalue itself. One can equally handle it numeri- 
cally by writing 



<5ct 4™ = (bn)^- ^ctHTJd' bri- 



(3.13) 



Rewriting i5ct in favor of renormalized parameters, 5^ — 
(5w/i;'™) (v™"(9/5w™" + V), yields 

Jeter = i^riSv/v'-n iVn+VTJn), 

dr 1 — 97' 

^ (1 + .^). 2.y a ^ (3.14, 

1 — ad 
for conciseness, we have suppressed "ren" in D. 
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FIG. 3: (a) Momentum profiles of the many-body corrections 
Ae'*~° /rj^'^"' for some typical channels, with v — l.lx lO^m/s, 
71 = 404 meV, A = 18meV and V4/V = -0.04 at B= 10 T. 



(b) 



-3l 



I Bo reveals the running of w"^"^" | b ■ 



With this in mind, let us rewrite Eq. 



+ (Ae''")'' 



(3.15) 



6a \rcn 



rcn 



where the renormaUzed correction (Ae ) 
(5ctea°" + Ae''° is now made finite. Writing the countert- 
erm as (^ctc™" - S^te 



if^" =T^b-Va+ Vir^b - Va), (3.16) 



and setting Ae''^°/77^'^" = Vc c''^", in units of the char- 
acteristic Coulomb energy 



V, = a/(efc£) « (56.1/efc) /bM meV, (3.17) 
yields the expression 

(Ae'"')"" = 77''^" {Vc c''^" + {V2/£) 6v}. (3.18) 
This reveals that the UV divergence is common to 



all ratios Ae"^"/??"^'' = Vcc"^", independent of 
(6, a) and (d, 1)4,71, A); the dimensionless quantities 
c''^'' have the structure c''^" = (V2/8) [log(AVB) + 
jT'ba^yren^ ^ren^ • • •)], where F''"' denote finite corrections. 

Figure 3 (a) shows for some typical resonance chan- 
nels the (rescaled) momentum profiles 7p [• • •] in Ae''*^° 
of Eq. p.7p . which, when integrated over i\p\, give 
^^b^aj^b^a ^^j^g p^^^g ^^Yai the slowly de- 

creasing high-momentum tail ^ [\/2/4)/{l |p|), common 
to all profiles, is responsible for the UV divergence. This 
numerically verifies the UV scaling law p.l8|) of the ra- 
tios Ae''^"/?/''^''. 

For renormalization let us refer to a specific channel 
(&o «o) and choose to define and other renormal- 
ized parameters by writing eoxcT"" = '^c"" {"Hbo ~ Vag ) at 
magnetic field B, or equivalently, (Ae^°"°)™" = 0, which 
yields Sv — — (a/\/2efc) c*""^"". The renormalized veloc- 
ity then runs with B, 



V2 



£6 



I So 



}, (3.19) 



and decreases gradually with increasing B. The leading 

-(V2/8) \og{B/Bo) 



correction c' 



bo-i-ao I 



is logarithmic but corrections coming from finite terms 
pboao g^j.g equally important for relatively low magnetic 
fields. For definiteness let us take L_2 <— L-3 as the 
reference channel, as chosen experimentally.^'^ For this 
channel the contribution from the low-momentum region 
decreases with B, as seen from the (—2 •(— —3) profiles 
for B = (10T,20T) in Fig. 3 (a), and numerically the 
correction is roughly doubled)^ 



-2< — 3 1 -2< — 3 1 
c Is - c \b„ 



^/2 , B 
-2.1 X — bg — 
8 ^Bo 



(3.20) 



over the range lOT < B < 30T, as shown in Fig. 3 (b). 
One can multiply it by factor a/{^/2ebv) ~ 1.5/eb ^ 0.3 
(with £{, ^ 5) to estimate the rate of decrease in v'^"^\b 
with B, which is about 10% for B = lOT ^- 20T. 

The renormalized Coulombic corrections in all other 
channels are thereby fixed uniquely. 



(Ae'"')"" = Vcrf^" {c' 



b<—a 



)\b. 



(3.21) 



These observable corrections are essentially calculated 
from the profiles in the low- momentum region £\p\ < 15. 
It is enlightening to write the resonance energies as 



^b-i—a 



Ac''''' = 



im - Va) 

„b-(—a 



VcAd 



b.a 



Vb - Va 



)U 



(3.22) 
(3.23) 



so that the Coulombic corrections Vc A c'''° seemingly 
arise relative to w"". Using the set of parameters in 
Eq. (|2.3p . one finds, for some typical intraband channels. 



Ac5'4 

A „3,2 1^=8 



Ac- 



0.317, (-0.259), 
= " -0.223, (-0.175), 
0.073, (-0.040), 



Ac^^i 0.650, (0.533), 

^^1,-2 -^^4 Qgg^^ (0.675), 

Ac-''-' -^-8 0, (0), 



A. 



,-4,-5 !^=^16 



-0.165, (-0.150), 
-0.275, (-0.247), (3.24) 



at B = lOT {B = 20T). For comparison, setting 
W4 = A = yields the electron-hole symmetric values 
Ac^'-' = 0.815, Ac-''-' = 0, Ac-''-^ = -0.159 and 
Ac-5 -4 = -0.303 at B = lOT. Similarly, some inter- 
band channels yield 

Ac''-' 0.353, (0.375), 

Ac''-' 0.370, (0.363), 

Ac''-^ ''==12 Q_261, (0.278). 



(3.25) 




1.3 
1.2 
1.1 

1. 
0.9 
0, 

1.4 



(c) 



i/=-l6 -4< — 5 
-3<--4 
v = -S -2<--3 



i/=-4 1< — 2 




FIG. 4: Cyclotron-resonance data of Ref. [23 with error bars, 
reorganized in the form ^atc^'/iv'' ~ V") ^-^d plotted in units 
of Wc = V2vo/£ (with vo = 1.1 x lO'^m/s) for 71 ^ 404 meV, 
A « 18meV and V4/V = —74/70 ~ —0.04. (a) Electron data; 
for clarity the data points, originally at B—{10, 12, 14, 16, 
18) T, are slightly shifted in B. (b) Theoretical expectation 
according to Eq. (|3.19p . with Vc/ijJc ~ 0.3 (or £5 ~ 5). (c) Hole 
data, (d) Theoretical curves. 



Note that those corrections are ordered regularly in mag- 
nitude for a sequence of resonances; see also the theoret- 
ical curves for e|;j^° in Figs. 4(b) and 4(d). 

Our formula p.22p summarizes the effect of renormal- 
ization in a concise form and is also useful in analyzing 
the experimental results. One may rescale the observed 
excitation energies e^^"" in the form eexc^/iVb — Va) and 
plot them in units of l/£ cx ^/B for each given value of B. 
The Coulombic many-body effect will then be seen as a 
variation in characteristic velocity [1 + 0(Vc)] from 

one resonance to another, and a deviation of w™" from 
^/B behavior would indicate the running of v'^™ with B. 



.5*-4 
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FIG. 5: Interband cyclotron resonances in bilayer inclusions 
in multilayer epitaxialgraphene reported by Orlita et al. Here 
some data from Ref. |2J are analyzed in the same way as in 
Fig. 4, using vo = 1.02 x lO^m/s and 71 « 385 meV for zero 
asymmetry A = 74 = 0. 



Figures 4(a) and 4(c) show such plots for the series of 
cyclotron resonance reported in Ref. 23. Both electron 
and hole data are plotted in units of Wc — \/2vQ/t 
40-\/b|tJ meV (with reference velocity — 1-1 x 10^ 
m/s kept fixed). These plots offer a closer look into the 
plots in Fig. 2. 

They are to be compared with Figs. 4(b) and 4(d), 
which illustrate how each resonance would behave with 
B, according to Eq. for Vc ~ ll^/Bff] meV (or 

Cb ~ 5). The = —8 curve for e^^^ ■^l 
ticular, represents the running of u'*^ 



|^=_8, m par- 
Is according to 
Eq. ((3TT9)) (normahzed to 1 at = 14T). These the- 
oretical curves and experimental data look similar but 
differ in details. They are not quite consistent, but there 
are some notable features: (i) The 2 1 and 1 4— —2 
resonances (the v = ±4 data) appear distinct from the 
rest, especially in their variation with B. In addition, the 
f = ±4 theoretical curves are separated from the rest by 
appreciable Coulombic gaps, but such a gap is not seen 
in the hole data. This would indicate, as noted in Sec. II, 
that the n — (0, 1) sector in bilayer graphene is signifi- 
cantly modified from the naive one we have supposed. 

(ii) The ly — (8, 12, 16) electron data and the 1/ = 
(—8,— 12,— 16) hole data show a general trend to de- 
crease with B, consistent with possible running of 
with B. Such (~ logarithmic) running of v'^'^'^\b is a direct 
consequence of renormalization and is thus the key sig- 
nature of the Coulomb interaction. In both electron and 
hole data appears to run in the same way at a rate 
somewhat faster than naively expected. Such enhanced 
running could in part be attributed to possible quantum 
screening'^* of the Coulomb interaction in graphene such 
that £(, is effectively largcr^^^ for lower B. 

Interband cyclotron resonance was recently observed 
by Orlita et al.^'^ in bilayer inclusions in multilayer epi- 
taxial graphene on the C-face of SiC. They identify some 
rt — 1 -s— —n [or rt -s— — (n — 1)] resonances with n > 3 
and obtain, via fitting, u w 1.02 x 10® m/s and 71 ~ 385 
meV, which are somewhat smaller than those for bilayer 
graphene. 

Some of their data are analyzed according to our for- 
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mula (13.221) in Fig. 5; there we have set 74 = A = 
since this experiment searched for no intraband reso- 
nances which would clarify a possible electron-hole asym- 
metry. The data appear to indicate slight running of 
with B, far slower than in the data in Fig. 4 on 
bilayer graphene. This suggests that the Coulomb inter- 
action could be significantly weaker (or more efficiently 
screened) in multilayered epitaxial graphene than in ex- 
foliated bilayer graphene. 



of the zero-mode and pseudo-zero-mode sector in bilayer 
graphene. 
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IV. SUMMARY AND DISCUSSION 



Appendix A: Analysis of divergence 



Experiment suggests that bilayer graphene has intrin- 
sic electron-hole asymmetry due to subleading intralayer 
and interlayer couplings. In this paper we have studied 
cyclotron resonance in bilayer graphene with such asym- 
metry taken into account. 

The set of asymmetry parameters, A « 18 meV and 
74/70 ~ 0-04 derived from independent measurements, 
entails a considerable modification of Landau levels in 
bilayer graphene and improves the theoretical fit to the 
data on cyclotron resonance between higher levels |n| > 2 
in both electron and hole bands. In contrast, the fit to 

the data on (2 •(— 1) and (1 i 2) resonances appears 

somewhat puzzling, and this suggests that the zero- and 
pseudo-zero-mode Landau levels n = (0, ±1) are further 
affected by some sources other than A and 74. It would 
be important to clarify, both theoretically and experi- 
mentally, the detailed structure of this special sector in 
bilayer graphene. 

The Coulombic many-body corrections to cyclotron 
resonance in graphene, unlike in standard quantum-Hall 
systems, are afflicted with UV divergences, and one has 
to carry out renormalization to extract genuine observ- 
able corrections. We have shown how to perform renor- 
malization for bilayer graphene under a magnetic field 
by first constructing necessary conterterms in free space. 
This renormalization program, formulated analytically, 
can equally be handled numerically in practical calcula- 
tions by use of the reduced matrix Hamiltonian in 
Eq. JM]) and counterterm ScS^^ in Eq. (|3A3l) . As a 
further illustration, we present the renormalization pro- 
gram for monolayer graphene with a possible valley gap 
in Appendix B. 

Equation (|3.22[) summarizes the effect of renormaliza- 
tion on cyclotron-resonance energies in a neat and con- 
cise form. This formula is also useful in analyzing the 
experimental data; it magnifies possible effects of the 
Coulombic corrections per channel and running of the 
renormalized velocity 1;''°" with B, as we have seen in 
Sec. III. In particular, the nearly logarithmic running of 
u''"^" is a direct consequence of renormalization, specific 
to graphene.— 1^ It is remarkable that such a renormal- 
ization effect is apparently seen in the data. 

More detailed measurements of cyclotron resonance, 
both intraband and interband ones, are highly desired to 
pin down the many-body effects as well as the structure 



In this appendix we examine the UV structure of the 
Coulomb exchange correction. Let us first look at in 
Eq. (|2.ip and construct the electron propagator in free 
space, 

{^{x)^Hx')) = {x\z/iidt ~ n+)\x') (Al) 
with \x) = |x) \t). We divide %+ into 2x2 blocks. 



n+ = 



ma^ VP + V4Q I 

vP + ViQ'^ -mcr3 + A + 71CT1 ' ^ 



where we have set m = u/2; P — +pcr_ with 

P = Px + iPy and a± = (cti ± i(T2)/2\ Q = diag.(p^p); 
the unit matrix 1, which, e.g., multiplies A in Eq. (jA2|) . 
will be suppressed in what follows. 

We go to the Fourier (p, uj) space and invert {uj — H+) 
in this 2x2 block form. A direct calculation yields 
= ,7V^.^/D with 

Nil = T + {{v^ +vl)-fi+2v4v{Lu-A)}PaiP 

+m (S-2Aw + A2-72)a3, 
N22 = T -2v{vA-Viji)p^ 

+ij (7? - A2) + A (S + 2m2) 

—TO S (T3 + {71 (oj^ — m^) + 2 V4V uj} (Ti, 
N12 = v{uj^ ~ v+v^ + m^) P 

+2m U4 (T3P + Vi'E.Q 

+ {v4 ji — V A) {uj + 1710-3) P 

+ {v-(i-ViA){uj + ma-i)Q, (A3) 

and N21 ^ {Ni2)\ where 

r = u{{uj-A.f-{v''+vl)v'-m^~-il} 

+{(w2 + i;2) A-2w4V7i}p', 
S = cj^ + - TO^; (A4) 

v± = V and PaiP = {p^Y + P^ <^-- The denomi- 
nator D — det {lj — %+) is cast in the form 

D = {c^2_Ac^-(v2+t,2)p2_^2|2 

-(71 UJ + 2'!;4up^)^ 

-TO2(4w2p2-72_ a2)^ (A5) 
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which, for u = 2m 0, leads to the band spectra {ei} 
in Eq. ^T^. 

The Coulomb exchange correction to 0{a/£) is writ- 
ten as the convolution '}2,y^v\^iS{\) + k) of Vk = 27rQ;/|k| 
and the instantaneous limit of the electron propagator 
(**^)|t'=t = /(dcj/27r)(**t)^ p = j5(p). In partic- 
ular, divergences arise from the portion of S'(p), that 
decreases like l/|p| or slower for p oo, and we shall 
focus on that portion. 

Integration over w, with the standard boundary con- 
dition, is readily carried out, yielding, e.g.. 



du) 1 . a + b + c + d 

2^ D{uj) ^* {a + c){b + c){a + d){b + dy 



(A6) 



where D{u!) — {uj — a){uj — b){uj + c){uj + d) (a, b,c,d > 
0) is short for D in Eq. (|A5p . When the integrand is 
uj/D, replace the numerator on the the right-hand side 
with ab — cd; for uj'^ /D replace it with —ab(c + d) — {a + 
b)cd. Actually, this structure combined with the band 
spectra l|2.4l) gives rise to a simple rule to handle the 
integration over lo: One can effectively replace 



LjJ 
, ,2 



(A 

— w+t;- p 



r7i)/2 

2 I 



(A7) 



in TVjfe in evaluating their large-p behavior. 

One may further note the following: (i) The denomina- 
tor (o-|-c)...(6-|-d) ~ 16t;^D_|_w_ (p^)^ for large p. (ii) For 
W4 — )■ and A — )■ 0, -D(w) becomes an even function of 
uj. As a result, odd powers of to in N necessarily lead 
to W471 or A, and are one power of |p| less than naively 
expected; e.g., ab — cd^ 2{vA — 71) |p|; corrections of 
the form m do not arise since D is even in m = u/2. 

With this in mind one can now retain only the portion 

^11 « T + {{v'^ + vl)-ii + 2v4v{uj-l\)}P(JiP, 

N22 « r - 2i; (w A - 'y47i)p^ -f (710;^ + 2i;4i;p^ w) (Ti, 

N12 « v{ij^ -v+v-p^)P, (A8) 

for further consideration. Note that, in view of Eq. (|A7p . 
the particular combination 5, despite its appearance, 
yields no divergent correction. As a result, terms (x in 
Nil and N22 lead to finite corrections. This fact has the 
important consequence that the interlayer bias u = 2m 
requires no infinite renormalization. 

The structure of the 0{a/£) Coulombic selfenergy cor- 
rection precisely reflects the structure of S'(p) after con- 
volution with Vp (X l/|p|. Let us therefore compare 
Eq. with 'H+ in Eq. (|X2|) . Applying first the 

rule (IA7p to N12 in Eq. (jASp yields the asymptotic form 



{^^^)i2\t=t'=P/{2\p\) 



(A9) 



This P/2 IpI term leads to a divergent correction of the 
form oc P logA^, which requires renormalization of ve- 
locity^ V. Actually this leading form P/2 \p\ is the same 
as the one obtained previously^^ for U4 = A = 0, and 
this implies that (the divergent part of) velocity renor- 
malization is unaffected by the electron-hole asymmetry 
cx V4, A. On the other hand, the portion cx Q in (^'5'^)i2 



has the asymptotic form cx Q/p^ which leads to no di- 
vergence and this means that W4 remains finite. 

Let us next consider or iVn in Eq. (K8\i . Its 

PdiP portion has the asymptotic structure PaiP/\p\'^ , 
which, though potentially singular, actually yields no 
divergent correction via symmetric p integration [since 
PaiP oc or (p^)^]. The remaining portion cx F is com- 
mon to A^ii and N22- Those common F terms, though 
leading to a divergent correction, simply shift the zero of 
energy and are of no physical relevance. One can now 
eliminate this F term from (^4'^)22|t=t' and determine 
its asymptotic form 



22 t=t' 



V ji ~ V4 A 



vA 



■ 0-1 



V4 71 



4w+u_|p| ' 2w+u_|p| 

.(AlO) 

This implies that both 71 and A undergo infinite 
renormalization. Evaluating the convolution integral 
/ d^kwk ^(p + k) with momentum cutoff A eventually 
leads to the countcrtcrms in Eqs. p.9p ^ (|3.Iip . 



Appendix B: monolayer graphene with a valley gap 

In this appendix we outline the renormalization pre- 
scription for monolayer graphene with a possible valley 
gap M . The effective Hamiltonian is written as 



'H+ = vp ■ (7 + M as 



(Bl) 



at one valley and acts on a two-component spinor of the 
form — (V'AjV'b)*- One can pass to another valley by 
setting M -M and ^ ^f' = {-ip'g^ip'^Y . 
Using the instantaneous propagator 



, ^ 1 vp-a + Ma, 

2 ^1,2 p2 + M2 



(B2) 



one can calculate the Coulombic selfenergy correction to 
0{a) and find divergences of the form 



Sv 



^(a/Sefe) log A2, SM = 2 (Af""/w"") Sv, (B3) 



which implies that the mass gap Af , as well as v, under- 
goes renormalization. 

Let us now pass to the B case and denote the ze- 
roth order Landau-level spectrum as e^°" = wjj"" 7y„ with 
w;;™ = \/2 and 



V7 



= sgn[n] VH+P^^^*^J^- 



(B4) 



Letting Set = Sv ^/Si;'™ + SM d/dM'"' act on e^f" then 
yields the counterterm (5cte5^'"- In particular, one finds 
that I?77„ in Eq. p.l6p is now replaced by 



(B5) 



With Eqs. (jB4p and (jBSp . velocity and mass renormaliza- 
tion for monolayer graphene in a magnetic field is carried 
out according to formula (I3.22p : we have checked numer- 
ically that this renormalization program works correctly. 



10 



^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature (London) 438, 197 (2005). 

^ Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 
(London) 438, 201 (2005). 

3 Y. Zhang, Z. Jiang, J.P. Small, M.S. Purewal, Y.-W. Tan, 
M. FazIoUahi, J.D. Chudow, J. A. Jaszczak, H. L. Stormer, 
and P. Kim, Phys. Rev. Lett. 96, 136806 (2006). 

* N. H. Shon and T. Ando, J. Phys. Soc. Jpn. 67, 2421 
(1998); 

Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 
^ N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. 

Rev. B 73, 125411 (2006). 
^ W. Kohn, Phys. Rev. 123, 1242 (1961). 
'' For many-body corrections in quantum Hall systems, see 

K. Asano and T. Ando, Phys. Rev. B 58, 1485 (1998). 

Kohn's theorem is also discussed for graphene by R. 

Roldan, J.-N, Fuchs, and M. O. Goerbig, Phys. Rev. B 

82, 205418 (2010). 

* D. S. L. Abergel and V. I. Fal'ko, Phys. Rev. B 75, 155430 

(2007) . 

® A. Iyengar, J. Wang, H. A. Fertig, and L. Brey, Phys. Rev. 

B 75, 125430 (2007). 
^° Yu. A. Bychkov and G. Martinez, Phys. Rev. B 77, 125417 

(2008) . 

S. Viola Kusminskiy, D. K. Campbell, and A. H. Castro 
Neto, Euro. Phys. Lett. 85, 58005 (2009). 
K. Shizuya, Phys. Rev. B 81, 075407 (2010). 
" K. S. Novoselov, E. McCann, S. V. Morozov, V. L Fal'ko, 
M. I. Katsnelson, U. Zeitler, D. Jiang, F. Schedin, and 
A. K. Geim, Nat. Phys. 2, 177 (2006). 
E. McCann and V. L Fal'ko, Phys. Rev. Lett. 96, 086805 
(2006). 

T. Ohta, A. Bostwick, T. Seyller, K. Horn, and E. Roten- 
berg. Science 313, 951 (2006). 
E. McCann, Phys. Rev. B 74, 161403(R) (2006). 
1^ E. V. Castro, K. S. Novoselov, S. V. Morozov, N. M. R. 
Peres, J. M. B. Lopes dos Santos, J. Nilsson, F. Guinea, 
A. K. Geim, and A. H. Castro Neto, Phys. Rev. Lett. 99, 
216802 (2007). 

J. B. Oostinga, H. B. Heersche, X. Liu, A. F. Morpurgo, 
and L. M. K. Vandersypen, Nature Mater. 7, 151 (2008). 
M. Mucha-Kruczyriski, E. McCann, and V. I. Fal'ko, Solid 
State Commun. 149, 1111 (2009). 
^° Z. Jiang, E. A. Henriksen, L. C. Tung, Y.-J. Wang, 
M. E. Schwartz, M. Y. Han, P. Kim, and H. L. Stormer, 
Phys. Rev. Lett. 98, 197403 (2007). 

R. S. Deacon, K.-C. Chuang, R. J. Nicholas, 
K. S. Novoselov, and A. K. Geim, Phys. Rev. B 76, 
081406(R) (2007). 

E. A. Henriksen, P. Cadden-Zimansky, Z. Jiang, Z. Q. Li, 
L.-C. Tung, M. E. Schwartz, M. Takita, Y.-J. Wang, 
P. Kim, and H. L. Stormer, Phys. Rev. Lett. 104, 067404 
(2010). 

E. A. Henriksen, Z. Jiang, L.-C. Tung, M. E. Schwartz, 



M. Takita, Y.-J. Wang, P. Kim, and H. L. Stormer, Phys. 
Rev. Lett. 100, 087403 (2008). 

M. Orlita, C. Faugeras, J. Borysiuk, J. M. Baranowski, W. 
Strupihski, M. Sprinkle, C. Berger, W. A. de Heer, D. M. 
Basko, G. Martinez, and M. Potemski, Phys. Rev. B 83, 
125302 (2011). 

L. M. Malard, J. Nilsson, D. C. EUas, J. C. Brant, 
F. Plentz, E. S. Alves, A. H. Castro Neto, and M. A. Pi- 
menta, Phys. Rev. B 76, 201401(R) (2007). 
L. M. Zhang, Z. Q. Li, D. N. Basov, and M. M. Fogler, Z. 
Hao, and M. C. Martin, Phys. Rev. B 78, 235408 (2008). 
Z. Q. Li, E. A. Henriksen, Z. Jiang, Z. Hao, M. C. Martin, 
P. Kim, H. L. Stormer, and D. N. Basov, Phys. Rev. Lett. 
102, 037403 (2009). 

A. B. Kuzmenko, E. van Heumen, D. van der Marel, P. 
Lerch, P. Blake, K. S. Novoselov, and A. K. Geim, Phys. 
Rev. B 79, 115441 (2009). 

J. Nilsson, A. H. Castro Neto, F. Guinea, and N. M. R. 
Peres, Phys. Rev. B 78, 045405 (2008). 
D. S. L. Abergel and T. Chakraborty, Phys. Rev. Lett. 
102, 056807 (2009). 

T. Misumi and K. Shizuya, Phys. Rev. B 77, 195423 
(2008); K. Shizuya, Phys. Rev. B 75, 245417 (2007). 
Numerically, for u = 5meV crossing e|„|=i ~ eo+ takes 
place at B « lOT. 

Y. Barlas, R. Cote, K. Nomura, and A. H. MacDonald, 

Phys. Rev. Lett. 101, 097601 (2008). 

K. Shizuya, Phys. Rev. B 79, 165402 (2009). 

R. Cote, J. Lambert, Y. Barlas, and A. H. MacDonald, 

Phys. Rev. B 82, 035445 (2010). 

S. M. Girvin, A. H. MacDonald, and P. M. Platzman, Phys. 
Rev. B 33, 2481 (1986). 

J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Nucl. 
Phys. B 424, 595 (1994). 

For an eigenmode |^„) of _ff with eigenvalue e„ one gener- 
ally finds that (-^nl^ffli/jn) = 5tn. for an arbitrary variation 
5. 

c''"'*""'' |s is calculated by integrating the underlying mo- 
mentum profile over the interval < ^ |p| < A/V e_B, with 
the momentum cutoff A kept fixed. 

Recently the effect of renormalization on the electronic 
spectrum has been measured for suspended graphene, in 
which v^°^^ runs with the electron density p as u'^'^"|p ~ 
■v"""\pa + (a/8e6) log(p/po). See, D. C. Elias, R. V. Gor- 
bachev, A. S. Mayorov, S. V. Morozov, A. A. Zhukov, 
P. Blake, K. S. Novoselov, A. K. Geim, and F. Guinea, 
preprint arXiv:1104.1396 

Renormalization of electron interactions in bilayer 
graphene under zero magnetic field was discussed by Y. 
Barlas and K. Yang, Phys. Rev. B 80, 161408(R) (2009); 
O. Vafek and K. Yang, Phys. Rev. B 81, 041401 (R) (2010); 
R. Nandkishore and L. Levitov, Phys. Rev. B 82, 115431 
(2010). 



